STABILITY OF TRAVELLING WAVES IN 
STOCHASTIC NAGUMO EQUATIONS 
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Abstract. Stability of travelling waves for the Nagumo equation 
on the whole line is proven using a new approach via functional in- 
equalities and an implicitely defined phase adaption. The approach 
can be carried over to obtain pathwise stability of travelling wave 
solutions in the case of the stochastic Nagumo equation as well. 
The noise term considered is of multiplicative type with variance 
proportional to the distance of the solution to the orbit of the 
travelling wave solutions. 



1. Introduction 

The purpose of this paper is to introduce a new approach to the study 
of (local) stability of travelling waves and pulses in excitable media 
that is in particular well-suited for stochastic perturbations. We are 
interested in the classical equations modelling the propagation of the 
action potential travelling along the axon of a neuron. As a starting 
point in this paper we consider the Nagumo equation on the real line (cf. 
[10J) perturbed by stochastic forcing terms. We make particular use of 
the explicit knowledge of the travelling waves in this case. However, our 
approach will be robust w.r.t. small perturbations in the coefficients. 

Since the spectral considerations, employed in the classical stability 
analysis of nerve axon equations (cf. [3j [5j |6] and the recent monograph 
[2]) are not easy to carry over to the stochastic case, we look for a path- 
wise stability analysis in the sense of the classical Lyapunov approach 
to the stability of dynamical systems. A first novelty of the paper is the 
introduction of an additional dynamics of gradient type that adapts a 
given solution of the stochastic Nagumo equation to the correct phase 
of the travelling wave. This explicitely given phase adaption, which is 
in addition easy to implement numerically, is the analogue of the phase 
conditions introduced as algebraic constraints in the classical stability 
analysis (see in particular |5J). As a second novelty in this paper, we 
replace the usual spectral considerations, applied to the Schrodinger 
operator, obtained as linearization of the underlying dynamics along a 
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given travelling wave, by functional inequalities of Poincare type. Our 
hope is that the latter method will be generalizable also to general 
systems of reaction diffusion type because it only uses partial infor- 
mation of the travelling wave solutions. Certainly, it is well suited for 
stochastic perturbations as demonstrated in this paper. An additional 
advantage is that, in contrast to the usual spectral considerations, our 
approach allows explicit quantitative estimates, both, in the determin- 
istic and in the stochastic case and sensitivity considerations w.r.t. the 
coefficients. 

The paper is organized as follows: In Section |2] we first present our new 
approach in the case of the deterministic Nagumo equation, to demon- 
strate the main arguments in a somewhat easier setting. The analogue 
to the usual spectral considerations of the Schrodinger operator asso- 
ciated with the linearization along a travelling wave is contained in 
Theorem 12.31 Theorem 12.61 then contains our main result on the lo- 
cal stability of travelling wave solutions. In Section |3] we consider the 
Nagumo equation perturbed with multiplicative noise. The variance of 
the noise term is chosen to be proportional to the distance of the solu- 
tion w.r.t. orbit generated by the travelling wave of the deterministic 
Nagumo equation. Hence, combining our stability analysis of Section 
[2] with a careful analysis of the stochastic perturbation, we obtain in 
Theorem 13.11 the stochastic analogue of our local stability result in the 
deterministic case. 

In addition to the new approach to the stability analysis via functional 
inequalities we also would like to mention that the type of stochastic 
Nagumo equations considered in this paper are also new in comparison 
with the models of spatially extended neurons subject to noise studied 
numerically and analytically by Tuckwell and Jost in [T2| [13] and also 
by Lord and Thiimmler in [9J . In order to ensure existence and unique- 
ness of a solution to our stochastic partial differential equation we use 
the variational approach to stochastic evolution equations as presented 
in monograph [TT] with recent extensions presented in [8]. In particular, 
we make use of the Ito formula, that can be obtained for the Hilbert 
space norm of the variational solution. The implied semimartingale 
decomposition can then be used to apply the one-dimensional (time- 
dependent) Ito formula to any smooth transformation of the Hilbert 
norm. 

2. The deterministic case 
Consider the Nagumo equation 

(1) d t v(t,x) = vd 2 xx v(t,x) + bf(v(t,x)) (t, x) G R+ x R 
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on the real line, with v, b > and 

f( v ) =v{l-v){v-a) ae (o,^ 

The equation is obtained from the well-known Fitz-Hugh Nagumo sys- 
tem 

d t v(t, x) = vd% x v(t, x) + bf(v(t, x)) - w(t, x) + I 
d t w(t,x) = e(v(t,x) -jw(t,x)) (t,i)el+xl 

by letting e \, 0, i.e., setting the recovery variable w constant, and fur- 
ther equal to the input current I. It is well-known that for parameters 
in the exitable region, the Fitz-Hugh Nagumo system admits a trav- 
elling pulse solution modelling signal propagation along the axon of a 
single neuron. The analogue for the Nagumo equation is a travelling 
wave front v(t, x) = v TW (x + ct), where 

(3) v TW (x)= (l+exp^-^* 

moving to — oo at constant speed c = a/2z/6 (| — a) (cf. pQ). We are 
interested in the local stability of this wave front in the function space 
H = L 2 (R). 

Before we can state a precise definition of stability, we need to in- 
troduce first our concept of a solution that we are working with. To 
simplify notations in the following we write v TW (t) = v TW (- + ct). Next 
(formally) decompose the function v(t, •) = u(t, ■) + v TW (t) w.r.t. the 
travelling wave. The resulting equation for u is then given by 
(4) 

dtuQt, x) = u%Ji(t, x) + b (f(u(t, x) + v TW (t)) - f(v TW (t))) 
(t,x) e R+ x R. 

For the precise definition of the Laplacian d 2 xx we need to introduce 
the Sobolev space V = H 1,2 (M) of order 1, equipped with the usual 
norm \\u\\y := J (d x u) 2 dx + \\u\\g. Clearly, V <-} H densely and 
continuously. Identifying H with its dual H' we obtain the embeddings 
V H = H' <-} V. Recall that w.r.t. this embedding the dualization 
v(fi u )v between / G V and u G V reduces to v'(f, u )v = (fi v )ii = 
J fudx, i.e. the scalar product in H in the case where / G H. The 
Laplacian d 2 xx then induces a linear continuous mapping A : V — > V, 
since y^Aw, v)y = — f d x ud x v dx < \\u\\v\\v\\y. 

The nonlinear term 

(5) G(t, u) := f(u(t, x) + v TW {t)) - f(v TW (t)) 

in equation (jlj) can be realized as a continuous mapping 

G( v ) : [0,oo] x V -)■ V 
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that is Lipschitz w.r.t. the second variable on bounded subsets of V 
with Lipschitz constant independent of t. Indeed, due to the elementary 
estimate ||w||oo — H M ||v> the Taylor representation 

G(t,u) = f(u + v TW (t)-f(v TW (t)) 

= f(v TW )(t)u + \f^\v TW (t))u 2 + l -f^{v TW {t))u* 
I o 

and uniform bounds on \\f^ k \v TW (t))\\ 00 , k = 1,2, 3, we have for w £ V 
that 

(G(t,u),w)< J \f(u + v TW (t))-f(v TW (t))\Hdx 

< ci||u||v (l + \\u\\%) \\w\W , 

hence 

(6) ||Gf(t,«)||^<ci|H| v (l + |H|y 
and 

(G(t, Ul )-G(t,u 2 ),w)< J \f( Ul + v TW (t))-f(u 2 + v TW (t))\\w\dx 

< C 2 (l + \\Ui\\y + ||« 2 ||v) ||«i - i^HirllHIv 

so that 

(7) \\G(t,ui) - G(t,u 2 )\\v < c 2 (l + ||wi||v + ||u 2 ||v) IK - u 2 \\h 

for finite constants C\ and c 2 depending on /itq 11 only. 

Note also that the sum i/Au + bG(t, u) satisfies the (global) monotonic- 
ity condition 

(8) (uAui + bG(t, Mi) — vAu 2 — bG(t, u 2 ), u\ — u 2 ) < brj\\ui — u 2 \\ 2 H 

where V = sup 5eM /'(0 = since (f(s) - f(t))(s -t)< V (s - tf 

for all s, t G R, and the coercivity condition 

(9) (i>Au + bG(t,u),u) < -u\\u\\v + (brj + v)\\u\\% 

since f(s)s = (f(s) - /(0))(s - 0) < r]s 2 for all s£K. 

It is now standard (see, e.g. Theorem 1.1 in [8J) to deduce for all u G H 
and all finite times T existence and uniqueness of a variational solution 
u G L°°([0, T] ; H) n L 2 ([0, T] ; V) satisfying the integral equation 

(10) u(t) = u + [ uAu(s,-) + bf(u(s)+v TW (s))-bf(v TW (s))ds 



associated with fll]). Clearly, we may consider this solution m as a 
solution on the whole time axes t > 0. 
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The integral J* vAu(s, ■) + bf(u(s)+v TW (s)) — bf(v TW (s)) ds appearing 
in the integral equation (TTUj) is well-defined as a Bochner integral in 
L 2 ([0,T] -V), since due to © 

' \\a u (s)\\ 2 vi + ||/(«( fl ) + ^(.)) - /(^(s))!! 2 ,, d s 
t 

2 J„ / 1 i _ Il„.f+Ml4 



< C / ||u(s)||frd5 1 + SUp < OO 

Jo \ te[o,T] J 

for all t > 0. In particular, the mapping £ i-)- u(t), [0, oo) — > V', is 
differentiable with differential 

^ = uAu{t) + &/(«(*) + v TW {t)) - f(v TW (t)) e V 
and therefore also locally Lipschitz. 

Definition 2.1. The travelling wave solution v TW is called locally stable 
in H if there exists 5 > such that for initial condition i>o with Vq — 
v TW G H and ||i>o — w Tiy ||_ff < 5 the (unique variational) solution u(t,x) 
to (jlj) satisfies 

hm \\ v (t,-)-v TW (t + t )\\ H = 
for some (phase) to £ 

The stability of travelling wave fronts for the Nagumo equation has 
been studied in many papers as a prototype example for metastability. 
The mathematical analysis of stability properties of v TW faces two 
major difficulties. The first one is the obvious fact that the reaction 
term f(u) in the equation ([1]) is not strictly dissipative in the sense 
that 

(vd 2 xx vi + 6/(«i) - vd 2 xx v 2 - bf(v 2 ),vt - v 2 ) < -k*\\ v i ~ v 2\\h 
for some > 0, or equivalently, the associated potential F(v) = 
fvo ftt) dt is n °t uniformly strictly convex, but a double- well poten- 
tial. This remains true if we fix v 2 to be equal to the travelling wave 
v TW or any of its spatial translates v TW (- — y). A first naive calcula- 
tion, exploiting the coercivity condition ([9]), only yields the following a 
priori estimate. 

Lemma 2.2. Let u E L°°([0, T); H) n L 2 ([0, T); V) be the unique solu- 
tion of (jlj). Then 

\\u{t)f H < e 2 ^\\u \\ 2 H \fte[0,T]. 
Proof. The coercivity condition (Q implies that 

j t \Ht)\\% = 2(uAu(t) + bG(t,u(t)),u(t)) 
<2br ] \\u(t)\\ 2 H te[0,T}. 
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Integrating up the last inequality w.r.t. t yields the desired inequality. 

□ 

However, restricting u to the orthogonal component of the derivative 
d x v TW of the travelling wave solution, i.e. j ud x v TW dx = 0, we have, 
under suitable assumptions on a the local dissipativity 

(udl x u,u) + b(f(v TW )u,u) < — K*\\u\\y + Ci\\u\\v\\u\\ 2 H 

for some > and some finite constant C\ which is implied by the 
following Theorem: 

Theorem 2.3. Assume that a > ~ ^1 — > hence K*(a) = y(l — 

2a) 2 < 1. Let 

= — : 7T— \ and C* = 3b\ — . 

6b + 40 1 + b(r] + 1) V 2v 

Then 

(ud 2 xx u + b(f'(v TW )u,u) < - Ksf \\ u \\ 2 v + a(u,d x v TW ) 2 
for all u G V . 

The proof of the Theorem is postponed to Section HJ 

The second difficulty in the mathematical analysis of the stability prop- 
erties of v TW is to identify the correct phase-shift of v TW (t + to) to 
which to compare the given solution v of ([I]) . To this end we introduce 
an auxiliary ordinary differential equation of gradient descent type as- 
sociated with the minimization of the distance between v and the set 
N = {v TW (• + C) | C e R} of all phase-shifted travelling waves. More 
precisely, given a solution v to the Nagumo equation ([[]) with initial 
condition vq satisfying vq — v TW G // 1,2 (1R), and given any relaxation 
rate m > (that will be specified later) we consider the ordinary dif- 
ferential equation 
(11) 

C{t) = -m(d x v TW (- + C{t) + ct), v TW (- + C(t) + ct) - v(t, -)) H 
C(0) = 0. 

The next Proposition states that the ordinary differential equation is 
well-posed. 

Proposition 2.4. Let v = u + v TW (t) be a solution to ([T]) with u G 
L°°([0,T\;H)nL 2 ([0,T\;V). Then 



B(t, C) = (d x v TW (- + C + ct),v TW (- + C + ct) - v{t, ■)) 



H 



is continuous in (t,C) G [0,T] x R, and Lipschitz continuous w.r.t. C 
with Lipschitz constant independent oft. 



Proof. Using the representation 



B(t, C) = (d x v TW (- + C + ct),v TW (- + C + ct)- v TW (t)) 
+ (d x v TW (- + C + ct),v TW (- + ct) - v(t, -)) H 



the continuity of B follows from the continuity of (t, C) >->■ d x v TW (• + 
C + ct) as a mapping with values in V, the continuity of (t, C) >->■ 
v TW {■ + C + ct) — v TW (t) as a mapping with values in V and the 
(Lipschitz) continuity of t v(t, ■) — v TW (t) as a mapping with values 



For the proof of the Lipschitz property w.r.t. C note that 



B(t,Ci)-5(t,C 2 ) 

= (a^ TH/ (- + d + ct) - d x v TW (- + C 2 + ct),v TW (- + ct) - v(t, •)) 
+ (d x v TW (- + d + ct),^- + d + ct) - v TH/ (- + ct)) H 
- (c^™(- + d + ct),^(- + d + ct) - t; TW (- + ct)) H 
= I + II + III, say. 



(13) 

\II + J//| = |(<9^ TW (- + ct),v TW (- + ct) - v TW (- - d + ct)) H 
- (cU™(- + ct),v TW (- + ct) - v TW (- - d + ct)>* | 

= + Ct),^(- - d + Ct) - - d + Ct)>n| 

< \\d x v TW (- + ct) \\ H 

■ fj II ex P(-\/J(- + t + c*))« 2W (- + e + ctf\\n d£ 




We next assume w.l.o.g. C\ < d- Then 
(12) 

^^(rc + d + ct) - v™> + d + ct)| 




so that we can further estimate 




< const • |d — d| • 



Similarly, using 
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fC 2 



r<^2 r 

\I\ < I I \d? x v TW (x + £ + ct)u(t,x)\ dxd£ 
(14) <3^£ 2 \\d x v TW (- + C + ct)\\ H \\u(t)\\ H d^ 

<hhr\Ci-C 2 \ sup \\u{t)\\ H . 



2u 



te[o,T] 



Inserting (lT5j) and (JHj) into ( II 2p yields the desired assertion. □ 

According to the last Proposition the function C defined by (TTT]) is 
well-defined. As already indicated, C will adapt to the correct phase 
of the v if we choose m > C* (cf. 12. 3p and our aim is to prove in the 
following that the difference 

(15) u(t) := u(t)+v TW (t)-v TW (-+C(t)+ct) = v(t)-v TW (-+C(t)+ct) 

converges to zero as t — > oo if the initial condition uq = vq — v TW 
is sufficiently small in the if-norm. In the next Proposition we first 
identify the resulting evolution equation for u. 

Proposition 2.5. Let u = v-v TW (t) G L°°([0, T]; i7)flL 2 ([0, Tj; V) be 
a solution of (jl]) and let u be defined by (|T5|) . Then u G L°°([0, T]; H) R 
L 2 ([0,T};V) again and u satisfies the evolution equation 

(16) _ 

nil — 

—{t) = uAu(t) + bG(t, u(t)) - C{t)d x v TW {- + C(t) + ct) 

= uAu(t) + bG(t,u(t)) 

- m(d x v TW (- + C(t) + ct),u(t))d x v TW (- + C(t) + ct) 

with 

G{t, u) = f{u + v TW (- + C(t) + ct)) - f(v TW (- + C(t) + ct)) . 
In particular, 

lf t \\m\%=-^\\dMmH+HG(t,m,m 

- m(d x v TW (- + C(t) + ct),u(t)} 2 . 

The proof of the Proposition is an immediate consequence of the prop- 
erties of v TW and the equations (J1J and (ITT]) . 

As usual we will now consider the linearization of the mapping G(t, u) 
around zero. To simplify notations, let v TW (t) := v TW (- + C(t) + ct). 
Then we can write 



(17) 



G(t,u) = f'(v TW (t))u + R(t,u) 



9 

where 

R(t, u) = f(u + v TW (t)) - f(v TW ) - f(v TW (t))u 

2 o 
satisfies the estimates 

(18) (R(t,u),u) < (A + a)\\u\\ 2 H \\u\\ v < (4 + a)||u||^||u||^ 
and 

(19) \\R(t,u)\\ v ,<(A + a)\\u\\ 2 H (l + \\u\\ v ) 

Similar to the classical stability analysis of the Nagumo equation we 
now use the information on the spectrum of the Schrodinger operator 
i , Au + bf'(v TW )u contained in Theorem 12. 31 with the above localization 
to obtain the first local stability result. 

2.1. Main result. 

Theorem 2.6. Assume that a > ~ ^1 — • Recall the definition 

of and C* in Theorem \2.3[ If the initial condition v q = Uo + v TW is 
close to v TW in the sense that 

\\u \\h <S- 



6(4 + a) 

for some 5 < 1 and v(t) = u(t) + v TW (t), where u(t) is the unique 
solution of (j3J), then 

|K*) - v TW (- + C(t) + ct)\\ H < e-^^Wvo - v TW \\ H . 
Here, C (t) is the solution of (fTTj) with m > C*. 

Proof. Let u(t) := v(t) — v (t) be as in ( TT5]) . Then Proposition 12.51 

and ( JT8|) imply that 

(20) 

\j t \\u{t) \\ 2 H = (vAu{t) + bf'(v TW (t))u(t),u(t)) + b(R(t, «(*)),«(*)> 

-m((d x v TW (t),u(t))) 2 
< {uAu(t) + bf'(v TW (t))u(t),u(t)} 
+ b(A + a)\\u(t)\\ H \\u(t)\\ 2 v - m {(d x v TW (t) , u(t))) 2 . 

Using translation invariance of uA and J {d x u) 2 dx. Theorem l2.3l vields 
the estimate 

{uAu(t) + bf'(v TW (t))u(t),u(t)} 

(21) ft n 2 

< -KJ\u{t)\\ 2 v + Cj u{t)d x v TW dx 
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Inserting §TQ) into (J2QD yields that 

H t Mt)\\ 2 H < -K*\\m\v + 6(4 + a)\\u(t)\\n\\u(t)\\ 2 v ■ 

In the next step we define the stopping time 

T:=m{*>0| K «)|| h >*^ 

with the usual convention inf = oo. Continuity of t i— > \\u(t)\\H 
implies that T > since ||lio||i? < ^bp^y- For t < T note that 

\j t Mt)\\ 2 H < "(1 " Q^Wmiv < "(I - S)n4u(t)\\ 2 H 
which implies that 

\\u(t)\\ 2 H <e-^-^\\u Q \\ 2 H 



for t < T. Suppose now that T < oo. Then continuity of 1 1- >■ \\u(1;)\\h 
implies on the one hand that ||w(T)||# = 8 b ^2+a) an< ^ on ^ ne °^ ner hand, 
using the last inequality, 



{T)\\ H = lim||w(t)|| H < e-^- 5 ^ T \\u \\ H < 8- 



t\r 6(4 + a) 

which is a contradiction. Consequently, T = oo and thus 

\\u(t)\\ H < e-^-^lluolln Vt>0 
which implies the assertion. □ 



3. Stochastic stability 

We now turn to the stochastic Nagumo equation 

(22) dv(t) = [vdl x v{t) + bf(v(t))] dt + dW^t) 

with a multiplicative noise term. W = {W(t))t>o is a cylindrical Wiener 
process with values in some separable real Hilbert space U, defined on 
some underlying filtered probability space (Q, J 7 , (J-"(t)) 4 >o, P). The 
theory of Wiener processes on Hilbert spaces and associated stochastic 
evolution equations can be found in the monograph The disper- 
sion coefficient S, describing the standard deviation of the noise term, 
formally is a function of the L 2 -distance inf^e^ \\v — v TW (- + C)||# of 
v to the set Af of travelling waves, i.e., E(u) = S(t> — v TW (- + C)) for 
all Cel. 

As in the deterministic case we will give the equation a rigorous mean- 
ing by decomposing v(t) = u(t) + v TW (t) w.r.t. the travelling wave. 
The stochastic evolution equation for u is then given by 



(23) du(t) = [uAu(t) + bG(t, u(t))) dt + £(«(*)) dW (t) 
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where the nonlinear term G is as in fl5]) and we will now make the 
following assumptions on S. E : H i— > L 2 (U,H) satisfies 

(24) = E(u + w TW - w TW (- + C)) for any CgI, E(0) = 
and £ is Lipschitz continuous 

(25) I|s(^)-^2)|IL M <^IK-^|||. 

for some constant cr 2 . Here, L/2(U, H) denotes the space of all Hilbert- 
Schmidt operators L : U h-> H with the Hilbert-Schmidt norm ||-^||^ 2 ([/^) 

Y2k ll-^ e fc|l# f° r one (hence any) complete orthonormal system (ek)k of 
U. 

We consider the equation (123|) w.r.t. the same triple V <^-> H = H' 
V as in Section |2j Due to the properties ©, ©, © and ©, we can 
now deduce from Theorem 1.1. in |8] for all finite T and all (deter- 
ministic) initial conditions uq E H the existence and uniqueness of a 
solution (u{t))tE[QX\ °f ( f23|) satisfying the moment estimate 

E\ sup ||u(t)||jj r + / ||u(t)||y<ft) < oo 
\te[o,T] Jo J 

which implies in particular that u e L°°([0,T]; H)nL 2 ([0,T];V) P-a.s. 
As a consequence we can apply Proposition 12.41 to a typical trajectory 
u(-)(oj) to obtain a unique solution C(-)(uj) of equation (|ll|) . It is also 
clear that the resulting stochastic process (C(t)) t >o is (J-t)t>cr adapted, 
since (u(t)) t >o is. We will assume as in the deterministic case that the 
relaxation rate m is sufficiently large, i.e., m > C#. 

Similar to the deterministic case we now define the stochastic process 

u(t) = u (t) + v TW (t) - v TW (- + C(t) + ct) = v(t) - v TW (t) 

which is {J-t)t>Q adapted too and satisfies the stochastic evolution equa- 
tion 

du{t) = \uAu(t) + bG(t, u(t)) - C(t)d x v TW (t)\ dt + E(u(t) dW{t) 
where 

G(t,u) = f(u + v TW (t))-f(v TW (t)) 
and the moment estimates 

El sup \\u(t)\\ 2 H + [ \\u(t)\\ydt) <oo. 

\*G[0,T] JO J 

Due to [II], Theorem 4.2.5, we have the Ito-formula 

I ~ ||2 



u\\ H (t) = \\u(0)\\b+ I 2(uAu(s) + G(s,u(s)) 

2 

L 2 (U,H) 







C(s)d x v TW (s),u(s)) + m^))\\l 2 (u,H) ds + M t 
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with 

M t = 2 [ (u(s),E(u(s))dW(s)) . 
Jo 

It follows from the above representation that is a (scalar- 

valued) continuous local semimartingale, in particular we have also 
the (one-dimensional) time-dependent Ito-formula 

ft 

12 \ _ / Q,„/„ \]~J„\\]2\ , na , / ||~/„\||2 



<p(t,\\u(t)\\jj)= / + 29^(8, ||«(s)||^)(i/A«(s) 

Jo 

+ 6G(s, - C^S^^s), «(a)> 

(26) +^(a,||fi(a)ll^)l|E(fi(«))ll! 9 ^ 

+ ^^,||M( S )|| 2 )2||E*(«( S ))«( S )|| 2 a d S 

+ / ^(s, 115(3)11^) dM s 



for any y? G C 1,2 ([0, T] x Here, E*(u) denotes the Hilbert-Schmidt 
operator from H to U adjoint to E(w). 



Theorem 3.1. Assume that a > | ^1 — yijjj ■ Recall the definition of 

andC* in Theorem \2.<ft and assume that a 2 < Letvo = Uq + v tw . 
Let v(t) = u(t) + v TW/ (£) ; where u(t) is the unique solution of the 
stochastic evolution equation (123"]) and tt(t) = w(£) + i> TW/ (£) — v TW {t). 
Let 

(27) T := inf{t > | \\u{t)\\ H > c.} , c. = - , 

with the usual convention inf = 00. Then 



p(t< 00) <- 2 \\m\\ 2 H 

ci 



and 

lim \\u(t)\\ H = 

t— >oo 

on the set {T = +00} 



Proof. Similar to the proof of Theorem 12.61 we have the following in- 
equality 

(pAu(t) + bG(t,u(t)) -C{t)d x v TW (t),u(t)) 
< -«,||T2(t)||^ + b{A + a)\\u{t)\\ H \\u{t)f v . 

In particular, 
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for t < T, where T is as in ( |27|) . Since also HE^t))]^^^ < 
c 2 ||"u(£)||h it follows from the choice of a that 
(28) 



Applying Ito's formula fl26|) to log(e + x), (128]) implies for t < T that 

12 



log(e+Ht)|||)<log( £ + ||S(0)|||)-f / " _ ds 

+ M 00_i( M < e >>* 

where 
It follows that 

£ + M < ( e + KOIIi) exp (-| jf + M« - i(M«)«) 

We can now take the limit e — > to obtain the pathwise inequality 
< llfi(O)llirejq) f-yt + M t - ^(M> 4 

where 

"* (u(s),£(ii(a))dW00) 



Hs)\\ 2 h 



Since T is a stopping time and exp(M t — |(M) t ) a nonnegative contin- 
uous local martingale, hence a supermartingale, the optional sampling 
theorem now implies that 



clE (e^ T l {T<oo} ) < limE (e^™\\u(T At)\\ 2 H 

< lim£ fexp(M T/U - J(M) TAt ) ) ||u(0)||| 

<\\m\\ 2 H- 

In particular 

P(T<oo)<l||S(0)||^. 
On the other hand note that on the set {T = oo} clearly 

|2 ^ii l|2 I K>± . - - I 



H*)Hh < l|fi(0)|||r exp ( -y* + M t - -(M) t J -> 

because (Mj)t>Q is a time-changed Brownian motion (see [7j, Theorem 
3.4.6 and Problem 3.4.7), i.e., there exists a Brownian motion (B t )t>o on 
a possibly extended probability space such that M t = B/m t , and thus 
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limsup^Mt - \{M) t = limsup t _ >00 (M) i (^) - \{M) t < +00, 
which implies the assertion. □ 

4. Proof of Theorem 12.31 

The proof of Theorem 12 . 3 1 requires a number of preliminary results. To 
simplify notations in the following we simply write v instead of v TW in 

the whole section. Let w(x) = v(l — v)(x) = ( 1H f e -L)2 with k = \J^- 
Proposition 4.1. Let u G C^(R) and write u = hw. Then 

(uAu + bf'(v)u,u) < — — J (d x h) 2 w 2 dx + — J h 2 w 2 dx 

< -(1 - «,(a))^ / (d x h) 2 w 2 dx + ^Qk(h,w 2 ) 2 . 



«„(a) = y(l-2a) 2 . 



Here 



Moreover, if a > max{~(l — y ^j) ; 0} ~ 0.226, then re* (a) < 1 and 
(z,A M + 6/>K M ><-^|^(l-re*(a)) ^ y (d x h) 2 w 2 dx+\\u\\ 2 H ^ 

Proof. First note that 

* pi / \ / A , „ Vw„, Vu> , \ 

//Am + bf'(v)u = uAh + 2v Vh + c /i iu 

\ w w ) 

because 

vAw + f'{v)w = cVw . 
Integrating against udx = hw dx yields 

(uAu + bf'(v)u,u) = j [vAh + 2z/— Vh \ hw 2 dx + c J h 2 Vwwdx 
= — v J {d x h) 2 w 2 dx — c J hVhw 2 dx 

V f ' \2„..2 , , C 2 /" . 2 2 



< — — J (d x h) w dx + — y /i iu dx . 
We will prove in Lemma [4.21 below that 



for k = y|- Using = «;*(a)|, we conclude that 



(uAu + bf'(v)u, u) < — (1 — «*(a))— ^ (d x h) 2 w 2 dx 

+ ~6fc hw 2 dx 

This proves the first inequality. 

For the proof of the second inequality observe that Lemma 14.21 again 
implies for any 7 G (0, 1) that 

a ...2 / / 1 -A^l „ /„\\ I/ i , : ) ...2 
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-(1-K*(a))- y (<9 x ./i) w < -(1 - 7)(1 - K*(a))- y (dj/i) rfx 

f (1 — kJo)) .. ||2 1/ (1 — K*(a)) 2 \2 

Choosing 7 = = 2Q 2 ^ 3fe yields the inequality 

<„A„ +bf(v)u,u)< -^-^MlL I (d,hf v/ dx + \\uf H 



Using the elementary estimate 



36(1 + 6\ 6fc < 36/JL 




66 + 40 4/ ~ V 2z/ 
we obtain the second estimate. □ 

The following Lemma has been used in the previous proof. 
Lemma 4.2. Let k := Then 

3 b 



(29) 



J h 2 w 2 dx < k J (d x h) 2 w 2 dx + Z 1 (^J hw 2 dx 

for all h G C*(K). Here, Z = j w 2 dx = is a normalizing constant. 
Proof. We will first show that 

(30) / {h - h(0)) 2 w 2 dx<K [ (d x h) 2 w 2 dx . 



To this end we will split up the estimate w.r.t. x > (resp. x < 0) 
and show that 

poo i*oo 

(31) / (h - h(0)) 2 w 2 dx < k (d x h) 2 w 2 dx 
Jo Jo 

and 

i-O i-O 

(32) / (h - h(0)) 2 w 2 dx < n (d x h) 2 w 2 dx. 
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Indeed note that for x > 0, using 

(h(x) - h(0)) 2 



J d x h(s)ds\ < J e ks ds J e- ks d x h 2 (s)ds 
' e kx -l) [ e- ks d x h 2 (s)ds. 



Integrating the last inequality against w 2 dx we obtain that 



/oo 1 poo poo —kx —2kx 

(h-h(0)) 2 w 2 dx<-j^ e- ks d x h 2 (s)j (1 + e _ fcx)4 dxds 



-L r e- ks d^h 2 (s) ( 1 - 1 3e , fc ' ] ,/.s- 



3fc 2 Jo " K " J V {l + e' ks ) 3 

10 1 f°° 



which gives flHT]) because k 2 = \- 

For the proof of (132|) note that w 2 (— x) = w 2 (x), so that (l3Tj) implies 
Clearly, combining (1ST]) and implies (15U|) . For the final step 
of the proof of the Lemma we define the probability measure /i(dx) := 
Z~ l w 2 dx. Then 

h 2 (x) n(dx) = Var ^(h) + I / /id/x 



< Z- 1 ^ (/i- /i(0)) 2 w 2 cfe + ^ hdfji 

< Z~ x k J (d x h) 2 w 2 dx + (^J hdfij 

which implies the assertion. □ 
Lemma 4.3. Let u G C^(1R) and write u = hw. Then 

Proof. Clearly, 

v i (d x u) 2 dx = v \ (d x h) 2 w 2 dx + 2v \ hd x hwd x w dx + v \ (d x w) 2 h 2 



— v J (d x h) 2 w 2 dx — v J h 2 d 2 x ww dx . 

Using the fact that 

—vd 2 xx w = bf'(v)w — cVw < (brj + ck) w 

and the obvious estimate ck = b (| — a) < b, we obtain the desired 
inequality. □ 
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Proof, (of Theorem 12. 3j) First let u G C^(M) and note that under the 
assumption a > | ^1 — yj^j, hence /c#(a) < 1, Proposition 14. II implies 



the estimate 

(uAu + bf'(v)u, u) < ~ q^^4q ( 1 ~~ K *( a )) f ^ y (9 x h) 2 w 2 dx + \\u\\ 2 H 

+ 3bJ^(u,d x v) 2 . 

Combining the last estimate with the previous Lemma 14. 3[ we obtain 
that 

(uAu + bf(v)u, < ~ 66 + 4Q ! + 6( * + i) ( y / ^") a ^+11*) 

+ 3bJ^(u,d x v) 2 . 



This implies Theorem 12.31 with 

ft* = — : 7T~~ ~; and C* = 36a / — . 

66 + 40 l + % + l) V 2u 

For general u E V we deduce the desired estimate by standard approx- 
imation of u with C^-functions in the K-norm. □ 
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